Granular fluid thermostatted by a bath of elastic hard spheres 
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The homogeneous steady state of a fluid of inelastic hard spheres immersed in a bath of elastic hard 
spheres kept at equilibrium is analyzed by means of the first Sonine approximation to the (spatially 
homogeneous) Enskog-Boltzmann equation. The temperature of the granular fluid relative to the 
bath temperature and the kurtosis of the granular distribution function are obtained as functions of 
the coefficient of restitution, the mass ratio, and a dimensionless parameter /3 measuring the cooling 
rate relative to the friction constant. Comparison with recent results obtained from an iterative 
numerical solution of the Enskog-Boltzmann equation [Biben et al., Physica A 310, 308 (202)] 
shows an excellent agreement. Several limiting cases are also considered. In particular, when the 
granular particles are much heavier than the bath particles (but have a comparable size and number 
density), it is shown that the bath acts as a white noise external driving. In the general case, the 
Sonine approximation predicts the lack of a steady state if the control parameter j3 is larger than a 
certain critical value /3 C that depends on the coefficient of restitution and the mass ratio. However, 
this phenomenon appears outside the expected domain of applicability of the approximation. 

PACS numbers: 45.70.-n, 05.20.Dd, 51.10.+y 
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I. INTRODUCTION 



The simplest model to describe the dynamics of gran- 
ular matter in the regime of rapid flow consists of an 
assembly of (smooth) inelastic hard spheres with a con- 
stant coefficient of normal restitution a [1]. The Liouville 
operator and the BBGKY hierarchy governing the time 
evolution of the phase space density and the reduced dis- 
tribution functions, respectively, can be extended to the 
case of dissipative collisions [2, 3]. By assuming that 
the pre-collision velocities of two particles at contact are 
uncorrelated (molecular chaos assumption), an approx- 
imate kinetic equation for the one-body velocity distri- 
bution function can be derived, thus extending the re- 
vised Enskog theory to the realm of dissipative dynamics 
[2, 4]. The Enskog equation accounts for spatial corre- 
lations through the equilibrium pair correlation function 
at contact g[n(r)] as a functional of the nonequilibrium 
density field n(r). On the other hand, in the special case 
of spatially uniform states the pair correlation function 
g(n) becomes a constant, so that the Enskog equation re- 
duces to the Boltzmann equation, except for an increase 
of the collision frequency proportional to g(n). 

The familiar concept of equilibrium is absent in a gran- 
ular fluid due to the collisional dissipation of kinetic en- 
ergy (which is transferred to the internal degrees of free- 
dom) by an amount proportional to 1 — a 2 . Even if the 
system remains in a uniform state (the so-called homo- 
geneous cooling state) [5] , the total kinetic energy mono- 
tonically decreases with time, unless some kind of ex- 
ternal forcing is acted upon the system to compensate 
for the collisional loss of energy and thus a steady state 



can be reached. From an experimental point of view, 
a fluidized steady state is usually achieved by violently 
shaking the container [6-15]. To mimic the effects of col- 
lisions with the vibrating walls in an otherwise uniform 
system it has become popular to assume that each parti- 
cle is subjected to a stochastic force with the properties 
of a Gaussian white noise [16-23]. The effect of this force 
is to produce frequent (and weak) random "kicks" to the 
particles between two successive collisions. On the other 
hand, the white noise force is not the only thermostatting 
mechanism proposed in the literature [24]. For instance, 
energy can be injected into the system by the action of an 
"anti-drag" force (Gaussian thermostat) [22, 25] or a con- 
stant force directed along the motion direction ( "gravity" 
thermostat) [22, 25]. The former thermostat is equiva- 
lent to a velocity rescaling in the freely cooling state [22] . 
On a different vein, Barrat et al. [26] have recently pro- 
posed a model in which the energy injected in vertically 
shaken granular systems is transferred to the horizontal 
degrees of freedom through collisions with an effective 
random coefficient of restitution. Regardless of the heat- 
ing mechanism, the common feature is that the granular 
fluid reaches a uniform nonequilibrium steady state char- 
acterized by a velocity distribution function /(v) differ- 
ent from the equilibrium Maxwell-Boltzmann (MB) dis- 
tribution 



/mb(v) = n ( 



/ m 



2ttT 
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(1.1) 



In this equation, m is the mass of a particle, n is the 
number density, and T is the granular temperature. The 
two latter quantities are defined in terms of the velocity 
distribution function as 



n = J dv/(v), 



(1.2) 
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where 

(„*) = I f dvv k f(v). (1.4) 

The deviation of / from /mb is usually monitored by the 
value of the fourth cumulant (or kurtosis) 

«-!$->• « 

A negative (positive) value of k implies that the distribu- 
tion is flatter (less flat) around the mean than the normal 
one and thus the distribution is said to be platykurtic 
(leptokurtic). In the cases of the white noise and Gaus- 
sian thermostats, the distribution is platykurtic (k < 0) 
for a > 0.7, while it is leptokurtic (k > 0) for smaller val- 
ues of the coefficient of restitution a [19, 22, 27]. On the 
other hand, the distribution is always platykurtic when 
the system is heated with the gravity thermostat [22]. 
As for the magnitude of k, it is typically smaller with 
the white noise thermostat than with the Gaussian and 
gravity thermostats. In the random coefficient of restitu- 
tion approach the value of k depends on the probability 
distribution of a (through the moments (a) and (a 4 )), 
but is otherwise positive definite [26] . Another important 
measure of departure from equilibrium is the high energy 
tail of the distribution. The asymptotic behavior of the 
distribution function for large velocities is generally of 
the form ln/(v) ~ — v a , where a = 1, a = |, and a = 2 
for the white noise, Gaussian, and gravity thermostats, 
respectively [19, 22, 27, 28]. In the case of the random co- 
efficient of restitution, the values of the exponent a seem 
to be highly dependent on the probability distribution of 
a, ranging from a w 0.8 to a rts 2 [26]. 

The lack of an equilibrium state in granular fluids is 
especially evident in the case of a mixture. According 
to the principle of equipartition of energy, a mixture of 
two fluids at equilibrium share the same temperature 
(Ti = T 2 ). On the other hand, a binary granular mix- 
ture in the homogeneous cooling state or in a noncqui- 
librium steady state exhibits two different granular tem- 
peratures [14, 15, 29-43]. In general, the temperature 
ratio T1/T2 depends on the mass ratio, the diameter ra- 
tio, the three coefficients of restitution involved, the mole 
fraction, and the volume fraction [29, 37]. Typically, 
the largest influence on T1/T2 is due to the mass ratio, 
the heaviest component having the largest temperature 
[14, 15, 37, 43]. An extreme example of this breakdown 
of the energy equipartition occurs in the homogeneous 
cooling of an infinitely heavy impurity particle in a sea 
of inelastic particles: if the cooling rate of the granular 
fluid is smaller than an effective impurity-fluid collision 
frequency, then the impurity/fluid temperature ratio di- 
verges [30, 32], an effect that can be characterized as a 
second-order nonequilibrium phase transition [32, 33]. 

In a recent paper, Biben et al. [44] have proposed an 
interesting alternative way of uniformly heating a gran- 
ular fluid to achieve a nonequilibrium steady state. The 



granular particles are assumed to be immersed in a bath 
of elastic particles kept at equilibrium at a certain tem- 
perature Tb. In the steady state, the energy loss due to 
inelastic collisions among the granular particles is bal- 
anced by the energy gain due to elastic collisions of the 
granular particles with the bath particles. The velocity 
distribution of the bath is given by the MB expression 
(1.1), except that n must be replaced by rib (number den- 
sity of the bath particles) and m must be replaced by nib 
(mass of a bath particle). This system cannot be consid- 
ered as a true binary mixture since the bath is assumed 
to be unaffected by the granular fluid. The role of the 
bath fluid is to provide a thermostatting mechanism for 
the granular fluid. However, this problem retains most 
of the basic features found in a mixture, such as the ef- 
fect of mass disparity and the competition between dif- 
ferent length and time scales. Although the bath is at 
equilibrium, the inelastic nature of the collisions among 
the granular particles prevent the steady state to be the 
equilibrium one [44]. First, the equipartition of energy 
is broken since the granular temperature T differs from 
the bath temperature Tb (actually, T < Tb). Second, 
the kurtosis k of the granular velocity distribution is not 
zero. Interestingly enough, the temperature ratio T /Tb 
and the kurtosis k are not only functions of the coeffi- 
cient of restitution a, but also of the mass ratio m/nib 
and a dimensionless parameter u> measuring the mean 
free path associated with the granular-bath collisions, 
relative to the one associated with the granular-granular 
collisions. In Ref. [44], Biben et al. solved numerically the 
(steady-state) Enskog-Boltzmann equation by an itera- 
tive method and obtained T/Tb and k for several choices 
of the parameters in the range < a < 1, m/rrib > 1, 
lu < 1. In addition, they observed that their numerical 
results were consistent with a Gaussian high energy tail, 
ln/(v) <~ — v 2 , although this asymptotic behavior was 
reached for extremely large velocities only [44]. 

The same system as in Ref. [44] has been studied by 
Barrat and Trizac [36] by approximating the granular 
velocity distribution function /(v) by its MB form (1.1). 
This allows one to obtain a closed cubic equation for the 
temperature ratio T/Tb- Comparison with the iterative 
numerical solution of Ref. [44] shows that the MB ap- 
proximation provides a good estimate of T/Tb for the 
values of m/rrib > 1 and lu < 1 considered [36]. On the 
other hand, the MB approximation is obviously unable 
to estimate the kurtosis k of the distribution function. 

The aim of this paper is to revisit the granular fluid 
thermostatted by a bath of elastic hard spheres pro- 
posed by Biben et al. [44]. The main goal is to ob- 
tain (approximate) algebraic expressions allowing one 
to get T/Tb and k in terms of a, m/mb, and u from 
the Enskog-Boltzmann equation in the (first) Sonine ap- 
proximation. This approximation consists of expand- 
ing /(v)//mb(v) — 1 in generalized Laguerre (or Sonine) 
polynomials and retaining the first term only, whose co- 
efficient is the kurtosis k. When this approximate form is 
inserted into the second- and fourth-order moment equa- 
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tions and terms nonlinear in k are neglected, one obtains 
a closed tenth-degree equation for the temperature ratio 
T/Tb and an explicit expression of n in terms of T/Tb . A 
similar method has been used by van Noije and Ernst [19] 
to estimate k as a function of a in the case of a granular 
fluid in the freely cooling regime or heated by a white 
noise thermostat, the results showing a good agreement 
with computer simulations [22, 27]. However, the situ- 
ation in the case where the thermostatting mechanism 
is provided by collisions with elastic particles is much 
more complicated. In particular, the a-dependence of 
k is strongly influenced by two independent control pa- 
rameters, namely the mass ratio m/m^ and the mean 
free path ratio to. 

The organization and main results of the paper are 
the following. The problem and notation are presented 
in Sec. II. To gain some insight into the different compet- 
ing scales of the problem, the MB approximation is also 
worked out in Sec. II. This suggests the use of the ratio 
(3 cx uj(1 + 77i/mb) 1//2 (l — a 2 ) between the cooling rate 
of the granular particles and the granular-bath collision 
rate (or friction constant) as a control parameter more 
convenient than the mean free path ratio u>. 

The Sonine approximation is carried out in Sec. III. 
As said before, this approximation yields a tenth-degree 
equation for T/Tb. A useful approximate solution to 
this equation is found by expanding the solution around 
the MB approximation and neglecting terms nonlinear 
in the deviation. In order to get more explicit results, 
some limiting behaviors arc considered, including the 
"colloidal" limit [m/nib — > oo, u> cx (m/mi,) -1 — > 0, 
(3 cx (m/mb) -1 / 2 — > 0] and the "white noise" limit 
[m/mt — ► oo, lu = finite, (3 cx (m/nib) 1 '' 2 — ► oo]. In 
the latter limit, one recovers the case of the white noise 
thermostat [19]. Section IV is devoted to a comparison 
with the numerical solution of Ref. [44] for < a < 1, 
m/rrib > 1, u> < 1. An excellent agreement is found, thus 
validating the reliability of the Sonine approximation, at 
least for lu < 1. 

The Fokker-Planck limit (m/my, — > oo) of the Enskog- 
Boltzmann equation is considered in Sec. V. An analysis 
of the high energy tail sho ws that for asymptotically large 
velocities (v ^> J 2Tb / m) the distribution function tends 
to a Gaussian, namely In / ss —raw 2 /2Tb. If the granu- 
lar temperature is much smaller than the bath temper- 
ature (T <C 7b), there exists an intermediate range of 
velocities y/2T/m <C v <C \J 2Tb / m where the distribu- 
tion function has the form of a stretched exponential, 
In / <~ — v 3//2 . In the white noise limit one has T/Tb — > 
and so the Gaussian form is pushed out to infinity and 
the stretched exponential becomes the only observable 
asymptotic behavior. 

Section VI shows that the Sonine approximation pre- 
dicts an interesting singular behavior. As (3 (or, cquiv- 
alcntly, lu) is increased for fixed values of a and m/my,, 
a critical value (3 c (a, m/m h ) is reached beyond which no 
steady-state solution exists. While for (3 < f3 c (a,m/my,) 
the temperature ratio T/Tb takes a well-defined sta- 



tionary value, it decreases with time (T/Tb - * 0) for 
(3 > (3 c (a,m/my,). Therefore, the Sonine approximation 
predicts the existence of a (first-order) phase transition 
when the friction constant becomes small enough as com- 
pared with the cooling rate and so the bath is unable to 
thermostat the granular fluid. Since this phenomenon 
takes place outside the domain of applicability of the So- 
nine approximation, it is not possible to assert in this 
context whether it is an artifact of the approximation or 
not. The paper ends with a summary and some conclud- 
ing remarks in Sec. VII. 



II. BASIC EQUATIONS. 
MAXWELL-BOLTZMANN APPROXIMATION 

Consider a fluid of inelastic, smooth hard spheres 
(species 1) of mass my = m, diameter o\ = a, and coeffi- 
cient of normal restitution an = a. The number density 
is n\ = n. The granular fluid is kept at a stationary, 
homogeneous, and isotropic state by a thermostat mod- 
eled as a bath made of a number density n 2 — ny, of 
hard spheres of mass m 2 = my, and diameter er 2 = cry, 
(species 2). The granular particles are assumed to col- 
lide elastically with the bath particles (i.e. a 12 = 1). In 
addition, the bath is assumed to be at equilibrium, un- 
affected by the granular fluid. Its velocity distribution 
function / 2 = /b is then 

AM-*(£gf"»(- ! 3r)- <"> 

Assuming the validity of the Enskog-Boltzmann de- 
scription, the velocity distribution function /1 = / of the 
granular fluid satisfies the equation 

dthiy) = Jii[vi|/i,/i] + Ji 2 [vi|/i,/ 2 ]. (2.2) 

The collision operators Jij[ v i\fi> fj] are [29] 

Jij[vi\fi, fj] = U.j "I j dv 2 J dBe(a- vi 2 )(S ; -vi2) 

x [^/iWl/itv^-^viJ/ih)], 

(2.3) 

where gij is the pair correlation function for particles of 
species i and j at the contact point r = <Tij = (ai+aj)/2. 
The primes denote pre-collisional velocities, 

v i = vi - Ha (1 + a' 1 ) (a ■ vi 2 )5, (2.4a) 

v 2 = v 2 + (l + a^ 1 ) (B ■ vi 2 )ct, (2.4b) 

where Vi 2 = vi — v 2 and fUj = mi/ (mi +rrij). The first 
term on the right-hand side of Eq. (2.2) is a nonlinear 
collision operator describing the rate of change of / due 
to the inelastic collisions among granular particles, while 
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the second term is a linear collision operator (Boltzmann - 
Lorcntz operator) accounting for the elastic collisions of 
granular particles against bath particles. Since mass is 
conserved in a single collision, we have the identity 



/ 



dvi Jij[\\\fi,fj] = 0. 



(2.5) 



On the other hand, energy is not conserved by collisions. 
Multiplying both sides of Eq. (2.2) by v\ and integrating 
over velocity we get 



d t (v 2 } =-(A 11 +A 12 ), 



(2.6) 



where 



An = ~f dv i vlJij[vi\fi, fj]. (2.7) 

The inelastic 1-1 collisions produce an energy loss, i.e. 
An > 0. As a matter of fact, the quantity An is pro- 
portional to the cooling rate of the granular fluid (: 



111 



3T/m 



(2.8) 



Collisions 1-2 are elastic, so the total kinetic energy of 
the two colliding partners is preserved. Therefore, the 
sign of A12 depends on whether T > Tb or T < T\>. Since 
in the steady state the mean kinetic energy of the gran- 
ular particles is smaller than that of the bath particles 
(T < Tb), the granular partner gains energy on average, 
so A12 < 0. In the steady state both terms exactly com- 
pensate each other and one has 



An 



Analogously, 



d t {v 4 



0. 



-(Bu+Bi 2 ), 



(2.9) 



(2.10) 



where Ba are the collisional moments 



Bij =~Jdvi vfJijivilfu fj}. (2.11) 
In the steady state, 

Bn+Bi2=0. (2.12) 

The main objective of this paper is to determine the 
granular temperature Ti—T and the fourth cumulant k, 
which are related to the second and fourth moments of 
the distribution function by Eqs. (1.3) and (1.5). To that 
end, I will use a first Sonine approximation in the next 
Section. Before that, it is instructive to get an estimate 
of T/Tb by means of the much simpler MB approxima- 
tion for /, as done in Ref. [36]. This also will allow us 
to introduce the two competing rates, along with some 
useful dimensionless parameters. If / is replaced by /mb 
[cf. Eq. (1.1)] in Eq. (2.7), one gets 



where [19, 29] 



A 



3 2T 



(o) 



(2.13) 
(2.14) 



4° 2 } =3 7 — x{\-x% 



2T 



In Eq. (2.14), 



2 {2T\ 1 ^ 2 
v = -V27 — na 2 gn (2.16) 

3 \m J 



(2.15) 



is an effective (self) collision frequency. In terms of it, 
the cooling rate (2.8) is approximated by Q C , where 



C (0) = v{l-a 2 ). 



In Eq. (2.15), 



1/2 



(2.17) 



(2.18) 



is a parameter measuring the temperature ratio, where 
M = M21 = (1 + m/m b ) _1 , and 



2 1/2 / 2T 
—n h o 12 gi 2 » - 



1/2 



(2.19) 



defines a characteristic rate for (cross) collisions of the 
granular particles with the bath particles. In the limit of 
one thermalized Brownian particle (n — > 0, m/m^ — ► oo, 
T = Tb), it is easy to show that d t (v) = — 7(v) , so 7 plays 
the role of a friction constant. Henceforth I will use this 
terminology when referring to 7. The cooling rate 
and the friction constant 7 define the two relevant time 
scales of the problem. Their ratio is 



p^ C — = 2-^ 2 u^ 1 / 2 {l-a 2 ), 
1 



where 



UJ 



na 2 gn 
nhO-j 2 gi2 



(2.20) 



(2.21) 



represents the mean free path associated with cross col- 
lisions, relative to the one associated with self collisions. 
The relevant dimensionless control parameters of the 
problem are the coefficient of restitution a, the mass ra- 
tio m/nib (or, equivalently, fi), and the dimensionless 
parameter ui, the latter encompassing the dependence 
on the densities and the diameters of the granular and 
bath particles. This is the parameter space considered in 
Ref. [44]. Alternatively, one can eliminate u> in favor of 
P and use a, m/m\, (or n), and (3 as control parameters. 
This will be the point of view adopted in this paper. 

The MB approximation for T/Tb consists of inserting 
Eq. (2.13) into the exact condition (2.9). This leads to a 
cubic equation for x, 



2x(x 2 -I) = 13, 



(2.22) 



whose physical root is 



■5 



| v^cos 



3 

^cosh 



■ sin 



■ cosh 



-1 (3^ 

1 



sin 



| sin" 1 (S#j)]} 



(2.23) 



The temperature ratio is then given, on account of Eq. 
(2.18), by 



T 



m/m b 



(1 + m/m b )x 2 — 1 ' 



(2.24) 



Thus the departure of x from 1 is a measure of the break- 
down of the cquipartition of energy. It is interesting to 
note that in this MB approximation the parameter x de- 
pends only on f3 oc lo[i~ x I 2 (\ — a 2 ), being independent of 
the two other parameters (a and m/m b ). As will be seen 
in Sec. Ill, this is not so in the Sonine approximation. 

Although Eq. (2.23) gives an analytical expression for 
the temperature ratio in the MB approximation, it is 
instructive to consider some limiting cases where Eq. 
(2.23) simplifies considerably. Let us consider first the 
limit where the cooling rate is much smaller than the 
friction constant, i.e. (3 — > 0. This limit includes as 
particular cases the quasi-elastic limit (a — ► 1 _ ), the 
limit of a vanishing mole fraction of granular particles 
(n/n b — > => lu — > 0), and the limit of small granular 
particles (cr/ob — > =^> ui — > 0). In cither case, Eq. (2.22) 
yields 



(2.25) 



and Eq. (2.24) becomes 



T 



l + -(m h /m + l)P 



l--(m b /m + l)0, (3^0 



P^O, (2.26a) 
(2.26b) 



Equation (2.26a) is more general than (2.26b), since 
the latter requires that m/m b is kept fixed, in which 
case the granular temperature is slightly below the bath 
temperature. However, if in addition to (3 — > we 
have m/m b -> 0, Eq. (2.26a) still holds, but (2.26b) 
only applies if the product (m b /m)(3 goes to 0. On the 
other hand, if a/a b — > and m/rrih oc (a/a b ) 3 , then 
(m b /m)(3 cx <T b /cr — > oo and Eq. (2.26a) yields 

T/T h w 2(m/TO b )/?' 1 oc er/cr b , cr/o- b -» 0. (2.27) 

In the opposite limit of a cooling rate much larger than 
the friction constant {[3 — > oo), one gets 



1/3 



/3 



oo 



and 



T 
71 



m/m b / 1 
1 + m/m b \2 



-2/3 



oo. 



(2.28) 



(2.29) 



A physically important hybrid limit corresponds to 
m/m b — > oo, er/(7b — * oo, and n/rib — > in the scaled 
way mjra b oc (cr/crb) 3 oc n b /n, so that the mass densi- 
ties of the individual granular and bath particles, along 
with their volume fractions are held fixed. This case was 
referred to as the "colloidal" limit in Ref. [44]. In this 
scaled limit wa/tsi m b /m and (3 oc ji 1 / 2 — ► 0, so that 
Eqs. (2.25) and (2.26b) apply again. 

As will be seen later, the MB approximation provides 
values for the temperature ratio in good agreement with 
those obtained from the iterative numerical solution of 
the Enskog-Boltzmann equation [44] . On the other hand, 
this approximation is unable to estimate the deviation of 
the distribution function from the MB form, as measured 
by the fourth cumulant (or kurtosis) k. This requires 
the use of a Sonine approximation, as done in the next 
Section. 



III. SONINE APPROXIMATION 

A. Description of the approximation 

The Sonine approximation [19, 29] worked out in this 
Section consists of expanding /(v)//mb(v) in Sonine 
polynomials and then truncating the series after the first 
non-zero term: 



/(v) w /mb(v) 



15 mv 
T ~ "2T" 



+ 



2T 



(3.1) 

This Sonine approximation is expected to be reliable as 
long as |k| remains relatively small (say |k| < 0.1). When 
the approximation (3.1) is used in Eqs. (2.7) and (2.11), 
A12 and B12 become linear functions of k, while An 
and Bu are quadratic functions of k. However, their 
quadratic terms arc further neglected, in consistency with 
the spirit of the Sonine approximation (3.1). Therefore, 



Ai i 



B (0) +B {1) n. 

13 v 



(3.2a) 



(3.2b) 



Inserting (3.2) into Eqs. (2.9) and (2.12), and after elim- 
inating k, one gets the following closed equation for the 
temperature ratio T/T b : 



A 



(0) 

11 



.(0) 
L 12 



B 



(i) 
11 



B 



B 



(o) 
11 



B 



A 



(i) 
11 



A 



(3.3) 
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Once T is determined, the cumulant k is simply given by 



K — — 



AO) 



A 



(0) 
12 



B 



(o) 
li 



B 



(0) 
12 



l(l) 
ill 



+ A {1) 



o(l) 
■°11 



+ -°12 



(3.4) 



Equation (3.4) provides an estimate of the cumulant k 
within the first Sonine approximation. Obviously, the 
smaller the value of the magnitude of the estimate (3.4) 
the more reliable it is expected to be. 

The expressions for A^{^ and B^' 1 ^ have been derived 
by van Noije and Ernst [19]: 



4(1) 



B (o) _ 
li — 



2T 



m 



16 11 ' 



(0) 



B n - 



2T 



m 



32 



(69 + 10a 2 ) + 



1-a 



.(0) 
in > 



(3.5) 



(3.6) 



(3.7) 



where A$ is given by Eq. (2.14). The quantities Af^ 
and B^' 1 ^ were evaluated by Garzo and Dufty in Appen- 
dices A and B of Ref. [29] from the Sonine approxima- 
tion in the general case of arbitrary coefficients of normal 
restitution an and aw Particularizing to an = 1 we 
get Eq. (2.15) and 



3 2T 



4 2 } = ^-M- 2 (4- 



3m) - tA 



(3.8) 



2TY 
m J 



x' 1 ^ - x 2 ) [8^ 4 + x 2 (5 - 8/i) + fi] , 

(3.9) 



= g"U~) ^ 5 [ 4x6 ( 10 - 33 M + 48/i 2 - 30m 3 ) 

+ fix 4 (64 - 87^ + 48/i 2 ) + x 2 fi 2 (9fi - 17) + 3^ 3 ] . 

(3.10) 

Given the values of the three independent parameters a, 
H = (1+m/mb) , and (3 (or, alternatively, a, /i, and u>), 
Eq. (3.3) yields an algebraic equation of tenth-degree for 
the unknown x. The high degree of this equation is due 
to the highly nonlinear dependence of the cross collisional 
integrals A12 and B12 on the temperature parameter x, 
as shown by Eqs. (2.15) and (3.8)-(3.10). The algebraic 
equation (3.3) has generally two or four negative real 
roots (which are unphysical) and two positive real roots, 
both larger than 1. The largest positive root diverges 
in the elastic limit (a — > 1), while the smallest real root 
tends to 1 in that limit. Therefore, the physical value 
of x is the smallest positive root of Eq. (3.3). Once x is 
determined, the temperature ratio is given by Eq. (2.24) 
and the cumulant k is obtained from Eq. (3.4): 



16a: 3 [2a;(a; 2 -l)-/3] 
3/3a; 3 + 4/i(4 - 3/z)x 2 - 4/z 2 



(3.11) 



The physical solution to the tenth-degree equation 
(3.3) is generally very close to the MB approximation 
(2.23). In particular, once ui is eliminated in favor of the 
parameter [3 defined in Eq. (2.20), the remaining depen- 
dence of x on a and \i is generally rather weak. On the 
other hand, it is precisely the difference between the MB 
and Sonine estimates of T/Tb, no matter how small it 
is, what is responsible for the prediction of a non-zero k. 
This is obvious from Eq. (3.11) and the MB approxima- 
tion (2.22). 

An approximate solution to Eq. (3.3) can be obtained 
by writing 



x w xq{\ + e), 



(3.12) 



where xq{[3) is the MB solution (2.23), and then neglect- 
ing terms nonlinear in e. The resulting expression for e is 
given in the Appendix. In consistency with the assump- 
tion (3.12), Eq. (3.11) becomes 



32^(3^-1) 



3/3x1 + Vl 4 - 3n)xl - 4^ 2 



(3.13) 



Equations (3.12) and 3.13), with xq given by Eq. (2.23) 
and e given by Eqs. (A.1)-(A.3), provide explicit expres- 
sions for the temperature parameter x and the kurtosis 
k as functions of a, /j, = (1 +m/mb) _1 , and (3. They are 
useful for practical purposes in those conditions where 
the MB approximation provides a good estimate of the 
temperature ratio and so n is relatively small. These 
include the physically relevant cases of m/rrih > 1 and 
u < 1 considered in Ref. [44]. When k is not small, the 
numerical solution of Eq. (3.3) can differ significantly 
from Eq. (3.12) and so the Sonine approximation may 
not be applicable any more. This point will be addressed 
more extensively in Sec. VI. For the sake of complete- 
ness, the results presented in the remainder of this paper 
correspond to the solution of the full Eq. (3.3), unless 
stated otherwise. 



B. Limiting behaviors 



In the limit (3 — > (i.e. a — > 1 
yields the asymptotic behavior 



or lj 



l + h/3, (3->0, 



where 



1 20 - (23 + 2a 2 V + 30^ 2 
16 5 - + 7[i 2 



0), Eq. (3.3) 
(3.14) 

(3.15) 



Although the coefficient h is not a constant, it is very 
close to the numerical value h — > j corresponding to 
the MB approximation (2.25). The temperature ratio, 
according to Eq. (2.24), is given by 



^ w [l + 2h(m h /m+l)/3] \ 
Jb 

w 1 - 2h{m h /m + 1)(3, (3 



(3^0, (3.16a) 
-> 0, (3.16b) 
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where in the last step it has been assumed that m^/m 
is finite. According to Eq. (3.11), the cumulant k is in 
this limit 



4/i-l 



(3, /?-(), 



(3.17) 



provided that the mass ratio is kept fixed. Equations 
(3.14) and (3.17) can be obtained as well from Eqs. (3.12) 
and (3.13), respectively, by making use of Eqs. (2.25) and 
(A.4). 

Equation (3.17) cannot be applied in the colloidal limit 
(3 oc fj, 1 ' 2 — > 0. In that case, it is necessary to evaluate x 
to second order in (3. The result is 



, 1 a 3 ( 39 



2a 2 ) 



B 2 



(3.18a) 



^i- V 3 ( 79 + 2a V 

T h 2 1 640 P 



The leading term in x coincides with that of the MB 
estimate (2.25). On the other hand, the sublcading term 
differs from — ^ and is needed to evaluate k from Eq. 
(3.11) with the result 



2a' 



20 



0- 



(3.19) 



Again, Eqs. (3.18a) and (3.19) can also be derived from 
Eqs. (3.12) and (3.13), complemented with Eqs. (2.25) 
and (A. 7). More in general, it can be easily checked from 
Eqs. (3.13) and (A.6) that Eq. (3.19) holds when (3 -» 
and it — > 0, regardless of the scaling law between (3 and fi. 
It is interesting to note that, according to Eq. (3.19), the 
velocity distribution is platykurtic for a > \/2/2 ~ 0.71, 
while it is leptokurtic for a < y/2/2. The same happens 
in the case of the homogeneous cooling state, as well as 
in the case of the white noise thermostat [19]. 

In the combined limit a /<7b — > with m/nib oc 
(^/fb) 3 , i-e. the granular particles are much smaller than 
the bath particles but the mass densities of the individual 
particles are comparable in both species, [or more gener- 
ally if (3 — > and m^/m — > 00 with (rrib/m)(3 — ► 00], Eq. 
(3.16a) yields 

T/Tb s=s (l/2/i)(m/rab)/3 _1 oc o/ob a/a h ^0, (3.20) 

where in that case fi — * 1 and so h = (27 — 2a 2 )/96. In 
this limit, Eq. (3.17) cannot be used and the cumulant is 
instead given by 



16(4/i-l) 8(3 -2a 2 ) 



5 + 2(4/i - 1) 63 - 2a 2 



ojob -» 0. (3.21) 



This value of k cannot be considered very small (Jj < 
K < f§)> so the Sonine approximation, while possibly still 
qualitatively correct, may not be quantitatively accurate 
in this limit. It is worth noting that Eq. (3.13) yields 
Kwf (4/i-l) in this limit /x -> 1, (3 -» 0. 



In the (non-colloidal) limit of heavy granular particles 
at constant u> (so n ~ rib, cr ~ o"b) we have m/nib — * 
oo^>^^0^/3^oo. In that case, the solution to Eq. 
(3.3) behaves asymptotically as 



where the amplitude A is 

16+ (45 + 6a 2 )(l - a) 



A 



'64+ (177 + 30a 2 )(l - a)' 



(3.22) 



(3.23) 



This parameter is very close to the numerical value A — > i 
given by the MB approximation (2.28). In this limit the 
temperature ratio T/Tb goes to zero as 



(3.24) 



(3 18b) anc ^ t ^ ie cummant tends to the (small) constant value 



k w H(2A-1) 



16(l-a)(l-2a 2 ) 
64+ (177 + 30a 2 )(l - a) 



, m/nib — » 00. 



(3.25) 



Interestingly, this expression coincides with the one de- 
rived by van Noije and Ernst for the white noise thermo- 
stat [19]. For this reason, I will call this limit the "white 
noise" limit. I will come back to this point in Sec. V. In 
this white noise limit, Eq. (3.13) becomes k s=s e/16, with 
e given by Eq. (A. 8). 

The expressions for T/Tb and k in the limits discussed 
above are summarized in Table I. 



C. Transient regime 

The Sonine approximation (3.1) not only allows one to 
get the steady-state values of the temperature ratio T/Tb 
and the cumulant k, but also their time evolution. Let 
us introduce the scaled time variable 



r(t) - / diV(i'), 
Jo 



(3.26) 



which counts time essentially as the average number of 
self collisions experienced by a granular particle. In that 
case, Eqs. (2.6) and (2.10) can be rewritten as 



d T T = — A*T, 
3 



d T n = -^-B* + t( K + i)A* 
15 3 



where 



_ Mi + A12 B * _ B11 + B12 



v(2T/m) 



v(2T/mf 



(3.27a) 



(3.27b) 



(3.28) 
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TABLE I: Summary of the expressions for the temperature ratio T/T b and the kurtosis k in the Sonine approximation for some 
limiting situations where /3 tends to or to oo. The parameter /3 is a measure of the cooling rate £ relative to the friction 
constant 7 and is a function of the coefficient of restitution a, the size ratio a/ah, the number densities n and n b , and the 
mass ratio m/m b [cf. Eqs. (2.20) and (2.21)]. The coefficients ft and A are given by Eqs. (3.15) and (3.23), respectively. Case 
1 represents the quasi-elastic limit. In cases 2 and 3 the mean free path associated with cross collisions is much smaller than 
the one associated with self collisions (w — » 0), but the masses of the granular and bath particles are comparable. In case 4 the 
granular particles are much smaller and much lighter than the bath particles. In case 5 the granular particles are much larger 
and much heavier than the bath particles but both species occupy a comparable volume fraction (colloidal limit). Finally, in 
case 6 the granular particles are much heavier than the bath particles but the mean free path ratio ui is finite (white noise 
limit). 



Case 


a 


a/a h 


n/n b 


m/mb 


P 


T/Tb 


K 

, ... /___ \'2 



1 ->• 1" fixed fixed fixed ~ (1 - a) 1 - 2h(m h /m + 1)0 (4ft -1) J 



2 fixed fixed ->0 fixed ~ n/n b -> 1 - 2ft(m b /m + 1)/3 (4ft - l) i*±s£aJ p 

3 fixed -^0 fixed fixed ~(a/a b f^0 1 - 2ft(m b /m + 1)/3 (4ft - 1) (1+m / m " )2 /3 



4 fixed -»0 fixed ~ {a/a b f -> ~ (a/a b ) 2 -> f a (m/m b )r' -> ^ 



5 fixed 00 - (a/a b )- 3 ~ (a/cr b ) 3 ^00 ~ (a/a b )-' i/2 ^ 1 - ±/3 ^f^/ 3 

6 fixed fixed fixed -^00 ~ (m/m b ) 1/2 -f 00 (A/3)^ 2/3 -f f (2A - 1) 



0.8 

b 




0.6 - 



K 




FIG. 1: Time evolution of the temperature ratio T/T b and 
the cumulant k, according to the Sonine approximation, for 
a = 0.5, w = 1, and m/m b = 1 (thick solid line), 10 
(dash line), and 100 (dash-dot line). The initial condition 
is T/T b = 1 and k — 0. The thin solid line in the upper 
panel represents the time evolution of the temperature ratio 
in the MB approximation for the case m/m b — 1. For the 
other two cases the MB curves are indistinguishable from the 
Sonine ones. 



Equations (3.27) are exact, but they do not form a closed 
set. On the other hand, if the approximation (3.2) is 
used, Eqs. (3.27) become a set of two nonlinear au- 
tonomous equations that can be solved numerically by 
specifying the initial values of T/Tb and k only. As an 
illustration, Fig. 1 shows the time evolution of T/T b and 
k for a — 0.5, u = 1, and m/rrih — 1, 10, and 100, start- 
ing from an initial condition with T/Tb = 1 and k = 0. 
We observe that after a few collisions per particle the 
steady state has been reached. Figure 1 also shows the 
time evolution of the temperature ratio in the MB ap- 
proximation for the case m/m^ = 1. For the other two 
cases the MB curves are indistinguishable from the So- 
nine ones, in consistency with the fact that k is rather 
small in those two cases. 



IV. COMPARISON WITH THE NUMERICAL 
SOLUTION OF BIBEN ET AL. [44] 

By an iterative numerical scheme, Biben et al. have 
solved the (steady-state) Enskog-Boltzmann equation 
(2.2) for several choices of the mean free path ratio w < 1, 
the mass ratio m/nib > 1, and the coefficient of restitu- 
tion a < 1. Figure 2 compares T/T b and k in the Sonine 
approximation with the numerical solution for u> = 1 
and m/mb = 1, 10, 100, and 1000. As can be seen, the 
agreement is excellent, except for a small tendency of the 
Sonine approximation to underestimate k for m/m^ = 1 
and large inelasticity. For the other three values of the 
mass ratio, the cumulant n is small enough to make the 
temperature ratio predicted by the MB approximation 
practically indistinguishable from that predicted by the 
Sonine approximation. The sequence of increasing values 
of m/mb at fixed uo in Fig. 2 leads to the white noise limit 
discussed in the previous Section. In fact, the asymptotic 
expression (3.25) describes very accurately the curve of 
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0.0 0.2 0.4 0.6 0.8 ^ 1.0 



a 

FIG. 2: Temperature ratio T/Tb and cumulant k versus a for 
u: = 1 and m/rrib = 1 (thick solid lines and circles), 10 (dash 
lines and squares), 100 (dash-dot lines and up triangles), and 
1000 (dot lines and down triangles). The lines correspond 
to the Sonine approximation, while the symbols correspond 
to the numerical solution of the Enskog-Boltzmann equation 
[44]. (The values of the temperature ratio for m/mb = 1000 
are not given in Ref. [44].) The thin solid line in the top panel 
represents the temperature ratio in the MB approximation 
for the case m/mb = 1. For the other two cases the MB 
curves are indistinguishable from the Sonine ones and are not 
plotted. 




0.0 0.2 0.4 0.6 0.8 1.0 



a 

FIG. 3: Temperature ratio T/T h and cumulant k versus a for 
(m/mb)u = 1 with m/m b = 10 (thick solid lines and circles) 
and 100 (dash lines and squares) . The lines correspond to the 
Sonine approximation, while the symbols correspond to the 
numerical solution of the Enskog-Boltzmann equation [44]. 
The temperature ratios predicted by the MB approximation 
are practically indistinguishable from the Sonine ones and are 
not plotted. The thin solid lines represent the asymptotic be- 
haviors given by the leading term of Eq. (3.18b) (top panel) 
and by Eq. (3.19) (bottom panel) for the case m/rrib = 10. 
For m/mb = 100 the asymptotic behaviors are indistinguish- 
able from the Sonine curves and are not plotted. 



K for m/mb = 1000. 

A sequence leading to the colloidal limit (u oc rrib/m, 
m/nib — > oo ) is shown in Fig. 3, where T/T, and k are 
plotted for fixed (m/m^LO = 1 with m/mb = 10 and 
100. Again, the agreement with the numerical solution 
is excellent. For m/mb = 10 and 100, the temperature 
ratio is already very well described by the asymptotic 
law T/T b = 1 - §/3 [cf. Eq. (3.18b)]. As for k, the 
asymptotic behavior (3.19) has been practically reached 
for m/mb = 100. 

Let us consider now the dependence of T/Tb and k 
on the mean free path ratio ui for fixed values of m/mb 
and a. Figure 4 shows the results for equal masses 
(m/mb = 1) and three values of the coefficient of resti- 
tution (a — 0.9, 0.5, and 0). Near the origin, (3 oc uu 
is small and then the asymptotic laws (3.16b) and (3.17) 
apply. As the inelasticity increases, the range of values of 
to for which the temperature decays linearly shrinks. On 



the other hand, since, as remarked by Biben et al. [44], 
the cumulant is practically a linear function of u, the 
asymptotic law (3.17) keeps being rather accurate even 
near u = 1. 

As is apparent from Figs. 2-4, the Sonine approxima- 
tion yields excellent estimates for the temperature ra- 
tio T/Tb and the kurtosis k for the cases considered in 
Ref. [44] . Even the MB approximation does a good job in 
estimating T/Tb, as first noted by Barrat and Trizac [36]. 
Because of this, the Sonine curves plotted in Figs. 2-4, 
which have actually been obtained from Eq. (3.11) and 
the numerical solution to Eq. (3.3), are practically indis- 
tinguishable from those obtained from the approximate 
solution (3.12) and (3.13). Does the good performance 
of the Sonine approximation mean that the distribution 
function is accurately described by the truncated expan- 
sion (3.1)? Not necessarily so if k is not very small, as 
illustrated in Fig. 5, where the reduced distribution func- 
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CO 



FIG. 4: Temperature ratio T/Tb and cumulant k versus u> 
for m/mb = 1 and (a) a — 0.9, (b) a — 0.5, and (c) a — 0. 
The thick solid lines correspond to the Sonine approxima- 
tion, while the symbols in the bottom panel correspond to the 
numerical solution of the Enskog-Boltzmann equation [44]. 
(The values of the temperature ratio for cases (a)-(c) are not 
given in Ref. [44].) The thin solid lines in the top panel rep- 
resent the temperature ratio in the MB approximation. The 
dotted lines in both panels are the asymptotic linear behaviors 
(3.16b) and (3.17). 



tion 



r(c)^n- 1 (2T/m) 3 / 2 /(v), 



r/y/Wfm, (4.1) 



is plotted for the case with the largest deviation from 
equilibrium considered in Ref. [44], namely equal masses 
(m/mb = 1), equal mean free paths (u> = 1), and to- 
tally inelastic collisions (a = 0) . Figure 5 shows that the 
Sonine approximation (3.1) partially accounts for the de- 
parture of the true distribution function from the equi- 
librium one. The remaining difference is due to Sonine 
polynomials of order higher than two. Notwithstanding 
this, the truncated expansion (3.1) is still useful for es- 
timating the collisional integrals (2.7) and (2.11). Of 
course, for other choices of the control parameters such 
that k is smaller than in the case considered in Fig. 5, 
the approximate distribution (3.1) is expected to be much 
more accurate. 



/ 




FIG. 5: Reduced velocity distribution function /*(c) for u> — 
1, m/mb = 1, and a = 0. The solid line is the numerical solu- 
tion of the Enskog-Boltzmann equation [44] , while the dashed 
line is the Sonine approximation (3.1) with k = 0.12984. The 



dotted line is the MB distribution /* (c) = tv 



-3/2 



FOKKER— PLANCK LIMIT. 

TAIL 



HIGH ENERGY 



A. Limit m/mb —> oo 

A physically important situation arises when the gran- 
ular particles are much heavier than the bath particles. 
This includes the colloidal limit [m/mb oo, w a 
m-b/m => (3 oc (mb/m) 1 / 2 — > 0], as well as the white noise 
limit [m/mb — ► oo, lu = finite => (3 oc (m/mb) 1 / 2 — * oo]. 
The predictions of the Sonine approximation in these lim- 
its have been analyzed in the previous Section. Nonethe- 
less, it is worth studying the limit m/nib — > oo at the 
level of the kinetic equation itself. 

If m/mb — > oo => fi w mb/m — ► 0, the Enskog- 
Boltzmann collision operator Ji2[v|/i, /y can be approx- 
imated by the Fokker-Planck operator (see, for instance, 
Appendix A of Ref. [30]), 



Ji2[v|/!,/ : 




v + 



m 9v 



/(v), (5.1) 



where 7 is given by Eq. (2.19). This confirms the in- 
terpretation of 7 as a friction constant. In the Fokker- 
Planck limit (5.1) the quantities A i2 and B12 can be eval- 
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uated exactly, 



2T T h T b 



(5.2a) 



B " = 15 -'(S) + < 5 - 2b > 

This means that in this limit the Sonine approximation 
I 



(3.2) for — (1,2) becomes exact. In fact, Eqs. (5.2) 
can be reobtained from Eqs. (2.15) and (3.8)-(3.10) by 
taking the limit fj, — ► 0, in which case x — ► {T^/T) 1 / 2 . 
On the other hand, the self collision operator Jn is ob- 
viously not affected by the limit m/m^ — > oo, so that 
Eqs. (3.2) for — (1, 1) are still approximate. In the 
MB approximation, the parameter x = (Tb/T) 1 / 2 keeps 
obeying the cubic equation (2.22). On the other hand, 
in the Sonine approximation the tenth-degree equation 

(3.3) simplifies to a quartic equation: 



320(l-a)x 2 (x 2 -l) + /? [64+ (177 + 30a 2 )(l -a)] x{x 2 -l)-lQQf3{l-a)x~2(3 2 [16 + (45 + 6a 2 )(l - a)] = 0. (5.3) 

I 



Moreover, Eq. (3.11) reduces to 



16 



— [2x(x 2 - 1) - 0\ . 



(5.4) 



The physical root of Eq. (5.3) can be obtained analyti- 
cally. Nevertheless, a simpler practical expression is pro- 
vided by Eq. (3.12) with e given by Eq. (A. 5). 

The above results apply in the limit m/mb — > oo re- 
gardless of the magnitude of (3. In particular, (3 remains 
finite if w — > as m/mb — > oo in the scaled way lo oc 
{m^/m) 1 / 2 . This may correspond to n^/n oc (<r/crb) 3 
and m/mb oc (cr/ob) 6 , so that the volume fractions of 
both species are comparable but the mass density of a 
granular particle is much larger than that of a bath par- 
ticle. Below the more interesting cases of the colloidal 
limit and the white noise limit (where (3 — > and (3 — > oo, 
respectively) are considered. 



B. Colloidal limit 

Let us assume that the heavy granular particles are 
made of a material with a mass density comparable to 
that of the bath particles, i.e. m/mb oc (er/ob) 3 - If, in 
addition, the partial volume fractions of both species arc 
comparable, we have n/n^ oc (a^/a) 3 . Consequently, 
u oc rrib/m and (3 oc (mb/m) 1 / 2 — ► 0. In that case, 
Eqs. (5.3) and (5.4) reduce to Eqs. (3.18a) and (3.19), 
respectively. Therefore, the same results are obtained in 
the colloidal limit following two routes: (i) first m/mb — > 
oo at fixed (3 and then (3 — > 0; (ii) m/mb — > oo and (3 — > 
simultaneously, with j3{m/m,h) 1 ^ 2 fixed. 



C. White noise limit 

The Fokker-Planck limit (5.1) assumes that m/mb — ► 
oo but otherwise the parameter u), essentially measuring 
the concentration ratio ro/rib, remains arbitrary. Now we 



In this situation where the concentration of the heavy 
granular particles is comparable to that of the bath par- 
ticles, the cooling rate of the granular fluid is much larger 
than the friction constant and, as a consequence, the 
steady state is reached with a granular temperature much 
smaller than the bath temperature. In this case, Eqs. 
(5.3) and (5.4) yield Eqs. (3.22) and (3.25), with A given 
by Eq. (3.23), and the temperature is given by Eq. (3.24). 
As said in Sec. Ill, the result (3.25) coincides exactly with 
that obtained in the case of a granular gas heated with 
a white noise thermostat [19]. To understand this, note 
that since T^/T >• 1, the drift term can be neglected 
versus the diffus ion term in the Fokker-Planck operator 
(5.1) (for v < yj2T h /m) with the result 



a 

2 



d_ 

dv 



/(v) 



(5.5) 



where 



e 



2T T h 
m \ 1 

= -C (0) A, 
m 



3/2 



(5.6) 



where in the last step use has been made of Eq. (3.24) 
and the definition (3 — C^/l- The right-hand side of 
(5.5) can be interpreted as representing the effect of an 
external white noise force F(t) with the properties [16- 
19] 



(F(t))=0, (F(t)F(f)> 



m 2 tf 



\8(t 



t% (5.7) 



consider that u> is finite, so that (3 oc (m/mb) 



1/2 



oo. 



where I is the 3x3 unit matrix. Now we see that this 
white noise thermostat can be obtained as a limiting case 
of a thermostat consisting of a bath of elastic particles. 
The stochastic force F(t) can be interpreted as arising 
from collisions of the granular particles with a compara- 
ble number of much lighter (and hotter) bath particles. 
The first line of Eq. (5.6) relates the strength of the force 
correlations with the size, mass, density, and temperature 
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of the bath particles, and with the size and mass of the 
granular particles [cf. Eq. (2.19)]. The second line holds 
in the steady state and relates £q with the properties of 
the granular fluid only. 



D. High energy tail 

The Fokker-Planck limit (5.1) allows us to speculate on 
the high energy tail of the velocity distribution function 
of the granular particles. Let us assume that 



/(v) « Ke~ 



v > y/2T/m, 



(5.8) 



where K, k, and a are unknown constants. If a < 2, it 
can be argued that [19, 28] 



JuW\f,f] ~ -Trna 2 gil vf, v > ^2T/m. (5.9) 

Inserting (5.8) into the Fokker-Planck operator (5.1), we 
have 



d ( T b 8 



m 

kv a > 1. 



(5.10) 



To fix ideas, let us consider those cases such that in the 
steady state T <C lb (i.e. (3 » 1). This happens, for 
instance, in the white noise limit. Then, there exists an 
intermediate regime of velocities with y/2T/m <C v <C 
^/2Tb/m. Let us assume (to be checked by consistency 
later) that in that case (Tb/m)akv a ~ 2 3> 1, so the bal- 
ance between Jn and J12 yields a = § and 



2 4irmna 2 gii 

K 



2V2tt 



r™f /2 _ (5.11) 



9T 7 (T b /T) 3 / 2 3A(l-a 2 )V2T/ 

where in the last step use has been made of Eqs. 
(2.16), (2.17), and (3.24). The condition kv a > 1 is 
consistently satisfied if v » y / 2T/m because kv a <~ 
(ttw 2 /2T) 3 / 4 . As for the consistency of the assumed con- 
dition {T\ } /m)akv a ^ 2 » 1, we have 



0-2 



m 



2\ 3/4 



2T 



mu \ 



(5.12) 



so this quantity is indeed large in the intermediate regime 
(mv 2 /2T) > 1 > (mv 2 /2T h ). In summary, we get 



/(v) w ife 



- Aw 



3/2 



\j2Tjm < w < y/2T h /m, (5.13) 



with fc given by Eq. (5.11). 

On the other hand, in the true regime of asymptotically 
large velocities, v » /3y / 27^/m, the term Jn can be ne- 
glected versus the diffusion and drift terms of the Fokker 
Planck operator J i2 and therefore a = 2, k = m/2T\>. 
Thus, 



/(v) « iTe""™ / 2Tb , w » [3^2Xj^i » v^T/^. 

(5.14) 



Strictly speaking, the asymptotic behavior (5.9) does not 
hold if a = 2, but it is still plausible that | Jn[v|/, /]| <C 
j^/T h /T(mv 2 /2T h )f if v > l3 v / 2T h /m (i.e. the most 
energetic granular particles are not affected by self col- 
lisions). The crossover from the stretched exponential 
(5.13) to the Gaussian decay (5.14) takes place around 
a velocity i>* that can be estimated from the condition 
mv 2 /2T\ ) <~ kvl' 2 , with fc given by Eq. (5.11). This yields 



2T 



1/2 



2V27T 



3A(1 - a 2 ) 



(5.15) 



This crossover velocity is indeed much larger than the 
thermal velocity because T <C 7b. 

The above analysis is not strictly valid if there is no 
a great disparity between both temperatures, so that 
the intermediate regime (5.13) does not exist. It is 
reasonable to conjecture that the high energy tail is of 
the form (5.8) with a = 2 but k ^ m/2Tb, since in 
that case Eq. (5.9) does not hold and one may expect 
that |Ju[v|/,/]| - v{mv 2 /2T)f instead. Moreover, if 
m > m b but the masses are not very disparate, the 
Fokker-Planck limit (5.1) is not valid cither. On the 
other hand, the results obtained here, along with the ev- 
idence of the numerical solution of Ref. [44] (see Fig. 8 of 
Ref. [44]) give support to the following scenario: For ve- 
locities v 3> y / 2Tb/m, the distribution function reaches 
an asymptotic form (5.8) with a = 2, i.e. it is a Gaus- 
sian; if, in addition, T < Tb, there exists a window of 
velocities much larger than the thermal velocity y / 2T/m 
but much smaller than y 7 2Tb / m where the asymptotic 
Gaussian form has not been reached yet and the distri- 
bution function can be approximated by the stretched 
exponential (5.13). 



VI. RANGE OF APPLICABILITY OF THE 
SONINE APPROXIMATION. CRITICAL 
BEHAVIOR 

It is not feasible to know a priori the range of validity 
of the Sonine approximation. In principle, it is expected 
to be valid as long as the value of |k| remains small. 
But how small must |k| be to trust the Sonine approx- 
imation? As comparison with numerical solutions have 
shown in Sec. IV, the Sonine approximation gives excel- 
lent results even if n ~ 0.1 (cf. Figs. 2 and 4). Since the 
Sonine approximation assumes that the velocity distri- 
bution function is close to the MB distribution, a con- 
sistency test is that the values for the temperature ra- 
tio obtained from both approximations are close each 
other. As a criterion to measure the range of applicabil- 
ity of the Sonine approximation, let us define a thresh- 
old value of the control parameter, (3^{a,m/m\,), such 
that if /3 < /3 t h (ct, m/m^), then the temperature ratio ob- 
tained in the Sonine approximation differs less than 10% 
from that obtained in the MB approximation. The values 
of the temperature and the cumulant at the threshold will 
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FIG. 6: The top panel shows the threshold value /3 t h versus 
a for m/rrih = 0.1 (solid line), 1 (dash line), and 10 (dot 
line). At a given mass ratio m/wib, the values of the temper- 
ature obtained from the Sonine and MB approximations differ 
less than 10% for those points lying below the corresponding 
curve. The middle and bottom panels show the temperature 
ratio and the cumulant, respectively, on the threshold curves. 



be denoted by T t ^(a, m/mb) and Kth.{ot, m/mb), respec- 
tively. Figure 6 shows the a-dependence of /3 t h, 7th/Ib, 
and K t h for m/mb = 0.1, 1, and 10. We observe that the 
threshold value /3th grows as the inelasticity decreases. 
Moreover, f3 t h is very sensitive to the mass ratio. As the 
granular particles become heavier, the range of values of 
the control parameter [3 for which the Sonine approxima- 
tion is expected to be reliable widens significantly. Figure 
6 also shows that the validity criterion (3 < /3 t h is com- 
patible with values relatively large for the cumulant n, 
especially for large m/mb- 
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FIG. 7: The top panel shows the critical value (3 C versus a for 
m/mb = 0.1 (solid line), 1 (dash line), and 10 (dot line). At 
a given mass ratio m/mb, the Sonine approximation predicts 
that there does not exist a steady state if (3 > /3 C - The mid- 
dle and bottom panels show the temperature ratio and the 
cumulant, respectively, on the critical curves. 



In the MB approximation, the asymptotic behavior 
(2.28) takes place when (3 — > oo, irrespective of whether 
this corresponds to m/mb — > oo or u — > oo. On the other 
hand, in the Sonine approximation the behavior (3.22) 
reflects the limit (3 — > oo provided that m/mb — > oo at 
constant u. In other words, Eq. (3.22) does not hold 
in the Sonine approximation at finite mass ratio when 
(3 — > oo (=> lo — > oo =>■ n/rib — > oo), i.e. in the limit of a 
vanishing concentration of bath particles. What happens 
then in that limit? 

The analysis of the tenth-degree equation (3.3) shows 
that, as (3 increases for fixed values of a < 1 and m/mb, 
a critical value /3 c (a, m/mb) is reached at which the two 
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positive real roots coalesce and no solution with posi- 
tive real x exists for (3 > (3 C . However, as (3 — > (3~ , 
the mathematical solution of the problem tends to well- 
defined values x c (a,m/mb) and n c (a,m/mb). Figure 
7 shows (3 C , T c /Tb, and k c versus a for the same val- 
ues of the mass ratio m/mb as in Fig. 6. We observe 
that at given m/mb, /3 C increases as the inelasticity de- 
creases (in fact, lim^!- /3 C = oo). Analogously, at 
given a, (3 C increases with the mass ratio m/mb- For 
(3 < {3 c (a,m/m h ), a linear stability analysis of Eqs. 
(3.27) shows that the steady-state solution (x, k) is lin- 
early stable. On the other hand, if (3 > (3 c (a, m/mb), 
a steady state is absent so k monotonically increases 
and T/Tb monotonically decreases with time. There- 
fore, the dynamical system (3.27) presents a (first-order) 
phase transition in the Sonine approximation: a stable 
fixed point exists for (3 < (3 c (a,m/mb) (with the "or- 
der" parameter T/Tb ^ 0), while no fixed point exists 
for (3 > (3 c (a,m/m h ) (with T/T h -> 0). 

It is difficult to assert whether the prediction of this 
phase transition is physically correct or it is just an arti- 
fact of the Sonine approximation employed. On the one 
hand, as the bath particles become more and more di- 
lute, so the friction constant becomes smaller and smaller 
versus the cooling rate, it is possible that a certain crit- 
ical value is reached beyond which the bath is unable to 
thermostat the granular fluid and, in addition, the veloc- 
ity distribution of the granular particles strongly deviates 
from the MB form [45] . On the other hand, the predicted 
critical behavior takes place clearly outside the range of 
validity of the Sonine approximation since the values of 
the cumulant at criticality, n c (a, m/mb), are too high, as 
seen in Fig. 7. These values are so high that the Sonine 
description cannot be quantitatively accurate. Whether 
it is qualitatively correct when predicting the existence 
of a critical behavior remains an open question. 

In any case, it is clear that a situation where the true 
k were small (so the Sonine approximation should then 
be reliable), while the Sonine prediction of n were large, 
is not self-consistent. Thus, one can assert that, irrespec- 
tive of whether a steady state with (3 > f3 c exists or not, 
the true kurtosis must be quite large, indicating that the 
granular distribution function must be highly distorted 
with respect to the MB distribution. 



VII. SUMMARY AND CONCLUSION 

The intrinsic nonequilibrium nature of granular fluids 
appears clearly in the case of uniformly heated systems. 
The system considered in this paper was originally pro- 
posed by Biben et al. [44] and consists of a fluid of in- 
elastic hard spheres (of mass to, diameter cr, and coeffi- 
cient of restitution a) immersed in a bath of elastic hard 
spheres (of mass mb and diameter Ob) kept at equilib- 
rium at a fixed temperature T\>. The granular particles 
are subjected to two competing effects: on the one hand, 
they experience dissipative collisions among themselves 



with an effective collision frequency v. on the other hand, 
they collide elastically against the bath particles with a 
characteristic rate (or friction constant) 7. The first ef- 
fect tends to produce a decrease of the granular parti- 
cle T with an associated cooling rate = v(\ — a 2 ), 
while the second effect tends to thermalize the granular 
fluid to the bath temperature 7b • When both mecha- 
nisms compensate each other, a steady state es reached 
with T < Tb (breakdown of energy equipartition) and a 
nonequilibrium velocity distribution with a non-zero kur- 
tosis k. The stationary values of T/Tb and k depend on 
the ratio of time scales (3 = /j, the mass ratio m/mb, 
and the coefficient of restitution a. 

In order to obtain T/Tb and n from the Enskog- 
Boltzmann equation, I have considered the (first) Sonine 
approximation, which is possibly the simplest approach 
allowing one to estimate the kurtosis of the distribution. 
As a consequence, k is expressed in terms of T/Tb, the 
latter quantity being given as the physical root of a tenth- 
degree equation. Explicit analytical results have been ob- 
tained in several limiting situations. In the limit of large 
friction constant relative to the cooling rate {(3 — ► 0) 
with finite mass ratio, as well as in the colloidal limit 
[m/mb — > 00, (3 cx (m/mb) -1 / 2 — ► 0], the breakdown 
of energy equipartition is weak (1 — T/Tb oc (3) and the 
distribution function is close to the Maxwell-Boltzmann 
(MB) form (k oc (3). However, if m/mb — > and (3 — > 
with (mb/m)(3 — > 00, as well as in the white noise limit 
[m/mb — > 00, /3 oc (m/mb) 1 / 2 — » 00], the steady-state 
granular temperature is much smaller than the bath tem- 
perature and the cumulant remains finite. 

In the limit of heavy granular particles (m/mb ^> 00) 
the Enskog-Boltzmann operator representing the colli- 
sions of the granular particles with the bath particles 
becomes a Fokker-Planck operator. Its high energy anal- 
ysis shows that the distribution function has a Gaus- 
sian tail In / w —mv 2 /2Tb for extremely large velocities 
(v > j3y/2T h /m). If (3 is large enough so that T < T b , 
th en ther e exists an intermediate window of velocities 
(y / 2T/m < v < y^2Tb/m) where the distribution func- 
tion has the stretched exponential form In / w —kv 3 / 2 , 
the coefficient k being given by Eq. (5.11). If the concen- 
tration and size of the granular particles arc comparable 
to those of the bath particles, then (3 cx (m/mb) 1 / 2 — ► 00 
and, consequently, T/Tb oc (m/mb) -1 / 3 — > 0. This im- 
plies that the drift term can be neglected versus the 
diffusion term in the Fokker Planck operator (for v <C 
y / 2Tb/m), so the effect of collisions with the bath par- 
ticles is indistinguishable from that produced by a white 
noise stochastic force. Thus the well-known case of a 
white noise thermostat is recovered from the more gen- 
eral case of a thermostat made of a bath of elastic hard 
spheres in the limit m/mb — » 00, while keeping fixed the 
concentrations and sizes of the particles. The fact that 
a granular fluid uniformly heated by a Gaussian white 
noise can be interpreted as a particular case of a gran- 
ular fluid heated by elastic collisions with particles of 
an equilibrium bath highlights the physical interest of 
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the latter system. While this system does not intend to 
represent faithfully the heating mechanism by vibrations 
usually employed in real experiments, it embodies most 
of the relevant nonequilibrium steady-state properties of 
a binary granular fluid. 

For finite mass ratio m/rrib and collision frequency ra- 
tio (3 one must solve numerically the tenth-degree equa- 
tion for T/Tb. Even so, this represents a formidable sim- 
plification of the problem with respect to the iterative nu- 
merical solution of the Enskog-Boltzmann equation car- 
ried out in Ref. [44]. Moreover, an approximate explicit 
solution to the tenth-degree equation can be obtained 
by neglecting nonlinear terms in the deviation from the 
MB solution. Comparison with the numerical data of 
Ref. [44] shows an excellent performance of the results 
obtained from the Sonine approximation. On the other 
hand, the numerical solutions considered in Ref. [44] were 
restricted to m/m b > 1 and (3 < 2~ 3 / 2 (l + m/m b ) 1//2 
and did not explore the whole parameter space. The re- 
sults presented in this paper show that, for fixed values 
of a and m/m bj the quantitative accuracy of the So- 
nine approximation worsens as (3 increases. For large /?, 
the distribution function is highly distorted from the MB 
distribution and so the kurtosis k is not small enough to 
make the Sonine method reliable. Quite interestingly, 
there exists a critical value (3 c (a,m/m,b) such that the 
Sonine approximation fails to provide a steady-state so- 
lution for (3 > (3 c {a,m/m\,). Since this phenomenon ap- 
pears beyond the domain of applicability of the Sonine 
approximation, it might be an artifact of the approxi- 



mation. Nevertheless, it is also possible that the Sonine 
prediction of a phase transition is qualitatively correct 
and the bath of elastic particles is unable to thermalizc 
the granular fluid if the ratio between the cooling rate 
and the friction constant is larger than a certain critical 
value [45]. The investigation of this possibility will be 
the subject of a future work. 
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APPENDIX: EXPRESSION FOR e 

If x = xq(1 + e), xq((3) being given by Eq. (2.23), is 
inserted into Eq. (3.3) and terms nonlinear in e are ne- 
glected, one gets a linear equation for e whose solution 
is 



where 



(1 - a) N(a,fj,,/3) 
32x 4 D(a,fi,P) 



(A.l) 



N(a, fi, (3) = [16/jlx% + (1 - 2a 2 - 16^)x 2 + 2/x] [64/z(l - fi)x 3 + 4/3(3 + 8^ - 6^ 2 )x^ + 12/3 2 x + 3/3 3 ] , (A.2) 

D(a,/j,,(3) = 240(l-a)^J 2 -6[2(l-a)(5a 2 + 52^) + 59 - 43a]4° 
+ [6(1 - a)(15a 2 + 223/i - 96^ 2 + 48m 3 ) + 227 - 99a] 

- [2(1 - a)^(343 + 46^ - 16/J 2 ) + 2(1 - a)a 2 (l5 + 16/x - 12m 2 ) + 33 - a] x% 
+2(1 - a) (j, [42 + 149^ - 34m 2 + 2a 2 (8 - 5m)] 4 

-2(1 - a)/i 2 (23 + 2a 2 + 18^)x 2 + 8(1 - a)fi 3 . (A.3) 



In the limit (3 — > with fixed a and m> one simply has 



h-^j/3, (A.4) 

where h is given by Eq. (3.15). On the other hand, if m — ► with fixed a and (3, we get 

_ 3^ (l-q)(l-2a 2 ) 

e ~ 8x1 10(1 - a) [16 - Za 2 {xl - 1)(1 - 2x1)} + ( x o ~ l ) I 127 + 354a; o - 3a ( 53 + 86a; o)] 

I 



(A.5) 



If both (3 and fi go to 0, Eq. (A.l) reduces to In particular, in the colloidal limit {(3 cx p}l 2 — > 0), one 

gets 

^^W + W (-A--6) 3(1 _ 2a2) 

£ 1280 - P ■ [AJ) 
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The same result is obtained by making x a — > 1 in Eq. Comparison with Eq. (3.25) shows that Eq. (A. 8) is close 
(A.5). toew (2A-l)/3. 

Finally, in the white noise limit (/? oc [i^ 1 ! 2 — > oo), Eq. 
(A.l) yields 

^ 1 (l-a)(l-2a 2 ) 

~ 4 16+(43 + 10a 2 )(l-a)' 1 ; 
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